Experimental measurement and Finite element modeling of extraordinary 
optoconductance in GaAs-In metal-semiconductor hybrid structures 
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We present a detailed discussion of extraordinary optoconductance (EOC). [Here optoconductance 
is synonymous with photoconductance.] Experimental data was acquired via macroscopic metal- 
semiconductor hybrid structures composed of GaAs and In and subjected to illumination from an 
Ar ion laser. A drift diffusion model using the finite element method (FEM) provided a reasonable 
fit to the data. EOC is explored as a function of laser position, bias current, laser power density, and 
temperature. The positional dependence of the voltage is accounted for by the Dember effect, with 
the model incorporating the excess hole distribution based on the carrier mobility, and thus the mean 
free path. The bias current is found to produce a linear voltage offset and does not influence the 
EOC. A linear relationship is found between the laser power density and the voltage in the bare and 
hybrid devices. This dependence is reproduced in the model by a generation rate parameter which 
is related to the power density. Incorporating the mobility and diffusion temperature dependence, 
the model directly parallels the temperature dependence of the EOC without the use of fitting 
parameters. 

PACS numbers: 72.20.Jv, 72.40.-fW, 72.80.Ey 



I. INTRODUCTION 



The resistance of any device can be divided into two 
components, one physical and one geometric. Contribu- 
tions to the physical component arise from factors such as 
the doping levels, impurities, bulk conductivity, and bulk 
mobility. The geometric contributions arise from factors 
such as the device dimensions, shape, and arrangement 
of the leads as well as any inhomogeneities that may be 
present. Most often, the physical contributions dominate 
the resistance, but by judicious choices, the geometric 
component can become the dominant component. The 
recently discovered EXX effectsi in metal-semiconductor 
hybrid structures (MSHs) of the type shown schemati- 
cally in Fig. 1(b) are primarily geometric in nature and 
rely on a reallocation of current between the metal and 
semiconductor regions due to an external perturbation. 
Here XX refers to magnetoresistance (MR), piezoconduc- 
tance (PC), etc. The perturbation can alter the interface 
between the metal and semiconductor regions or, as we 
show here, the bulk transport characteristics of either 
constituent. If the metal and semiconductor conductivi- 
ties differ by several orders of magnitude, as is the case 
with GaAs and In, small perturbations can dramatically 
change the flow of current through the device. As a re- 
sult, under perturbation, the measured voltage difference 
(depicted as VmSH in Fig. 1(b)) increases as compared 
to that of the metal or semiconductor alone. The large, 
often linear, change in voltage makes EXX devices an 
excellent choice for high resolution spatial sensing. Ad- 
ditionally, because they are non-magnetic, EXX devices 
can be used in applications where typical magnetic sen- 



sors are not suitablei^ 

The original EXX effect to be discovered was EMRi^ 
An InSb-Au MSH structure of macroscopic (and later 
nanometer) dimensions was fabricated using novel tech- 
niques described elsewhere^ The perturbation in this 
case was a magnetic field, redirecting the current through 
the semiconductor, whereas that current flows primarily 
through the metal with no applied field. A room temper- 
ature EMR of order 10^% was observed at fields of 5T. 
Theoretical modeling of EMR, using the finite element 
method (FEM) incorporating the effects of the bound- 
ary conditions at the metal-semiconductor interface, has 
shown excellent agreement with the experimental datai^ 
EMR has potential applications in a number of diverse ar- 
eas including ultra-high-density magnetic recording, po- 
sition sensing, and medical instrumentation. Scanning 
EMR probes have also been proposed as a way to con- 
currently measure topographic and magnetic field infor- 
mation with high spatial and temporal resolution^ 

The second demonstrated EXX effect was EPC.^ An 
InSb-Au MSH structure, similar to that used for EMR, 
was constructed. In EPC the perturbation is a uniax- 
ial tensile strain of the metal-semiconductor interface. 
Room temperature extraordinary piezoconductance as 
high as 500% has been demonstrated in the above de- 
scribed MSH structures Theoretical models based on 
both the FEM solution and the analytic solutions of 
Laplaces equation have accurately reproduced the EPC 
effect.^ EPC devices have obvious applications in the con- 
struction industry and have also been proposed for study- 
ing basic material properties (interface dynamics). 

Following the discovery of EMR and EPC it was real- 
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ized that an optical equivalent of the EXX phenomena, 
e.g. extraordinary optoconductance or EOC, should ex- 
ist, with photons providing the perturbation. Such a 
device would have applicability in the solar energy field, 
modern fiber optic relays, position sensors, and other op- 
toelectronic devices. We have recently provided a proof 
of principle demonstration of the EOC effect in GaAs-In 
MSH structures in a brief preliminary report >2i Here we 
give a full account of this effect and show a detailed FEM 
calculation that quantitatively accounts for the observed 
experimental results. 

A diagram of an MSH device is given in Fig. 1(b) 
with the connections shown schematically. The mesh 
shown in Fig. 1(c) contains fewer elements than the ac- 
tual mesh used for calculations. The bare device has the 
same physical connections as shown in Fig. 1(a). In a 
voltage measurement, the surface of the semiconductor 
region is exposed to a focused laser beam, whose lateral 
position (x£, yi) can be accurately controlled (see exper- 
imental details below). We define EOC as the percent 
difference in the measured output voltage in the MSH, 
as compared to that of the bare semiconductor with no 
shunt attached, such that 



[V23{xe, ye)]MSH 



EOC{x,,y,,X,{(3}) = 

ba', 



, X 100%. (1) 

[V23(Xe, y£j\bare 

Here the parameter set represents the power density 
of the laser, P, the temperature, T, and the bias current 
/. The quantity A = 476.5 nm is the emission wavelength 
of the laser, which is held fixed for all experiments. In 
this paper, we investigate the EOC as a function of the 
laser spot position with the parameters {/3} held fixed, 
and as a function of each of the {/?} parameters with 
the spot position and the other parameters in the set 
held fixed. Because of the many factors influencing EOC, 
it presents a unique opportunity to study the dynamics 
of the metal semiconductor interface, carrier transport 
phenomena, and the effect of introducing a shunt. 



II. EXPERIMENT 
A. Experimental Setup 

The experimental setup used to measure the EOC is 
shown in Figure 2. The bare (shuntless) sample was il- 
luminated with a Coherent Innova 400 Ar-t- ion laser op- 
erating at 476.5 nm with the beam focused to a 20 /im 
diameter. The power density was varied from 6.3 x 10"* 
W/cm^ to 5.8 X 10^ W/cm^. The samples were mounted 
in an ARS closed cycle helium cryostat and cooled to 
temperatures ranging from 10 K - 300 K. The cryostat 
temperature was maintained by a Scientific Instruments 
9650 temperature controller. The sample was positioned 
using a Newport Universal Motion Controller ESP 300 
system with three linear DC stepper motors. For op- 
timal spatial resolution, the sample was placed at the 



(a) 
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FIG. 1: (Color online) Bare (a) and MSHs (b) devices, (c) 
The FEM mesh showing the device dimensions, port configu- 
ration, and schematic of the connections to the device. 



focused beam waist, thereby defining the distance from 
the lens to the sample and the diameter of the laser spot. 
The bias current was supplied by a Lakeshore 120 cur- 
rent source in the forward and reverse direction from 1 
/J, A to 100 mA across ports 1 and 4, as shown in Fig. 
1. The voltage was acquired by a Keithley 2182 nano- 
voltmeter. One channel measured and the other channel 
(not shown) measured the voltage drop across the 1.2 Q 
resister in series with the sample. The voltage across the 
resistor was used to precisely determine the bias current 
flowing through the device. 



B. Sample Characterization 

Bare samples of dimension 2 mm x 10 mm x 0.4 mm 
were prepared by dicing a 2-inch diameter n-type GaAs 
wafer with the [001] growth direction along the 10 mm 
side. Gallium arsenide was chosen because it is a direct 
gg^piOji^ (^x.424 eV at 300 K) semiconductor and because 
of optimal absorption in the spectral region of the argon 
ion laser. The sample was degenerately doped"'^^ with Si 
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FIG. 2: (Color online) A schematic diagram of the experi- 
mental setup for sample positioning and data acquisition. 

at a concentration No = 1.25 x lO^^cm""^. The elec- 
tron mobility of the semiconductor was measured to be 
2100 cm^/Vs. Leads were attached first by metalizing 
the GaAs, in an inert nitrogen atmosphere, with dots of 
In (of diameter ^ 0.2 mm) arranged along the 10 mm 
side (see Fig. 2). Indium was chosen for metalizing be- 
cause of its low melting temperature and compatibility 
with GaAs.^^ This setup is similar to the typical four- 
lead van der Pauw plate setup used elsewhereii^ This 
geometry is the conformal equivalent^ to a van der Pauw 
disk with an off-centered inclusion. Next, the wires were 
tinned with In and pressed onto the metalized In dots. 
Ohmic contacts over the temperature range of interest 
were confirmed by verifying the linearity of the measured 
I-V response. These same bare samples were then used 
to make hybrid samples via the addition of an In shunt. 
This was achieved by metalizing the bare samples on the 
side opposite the leads with Indium of dimensions 10 mm 
X 0.5 mm x 0.5 mm. The hybrid sample resistance at 
room temperature, as measured across ports 1 and 4, de- 
creased by an order of magnitude from that of the bare 
samples, indicating the effectiveness of the shunt. Ohmic 
contact at the interface was again verified by the linearity 
of the I-V response. 



C. Experimental Procedure 

Once the sample position was calibrated, the positional 
dependence of the voltage V23 was studied by moving 
the laser spot along the x-direction for fixed values of 
yi. However, in no case was the interface between the 
semiconductor and the shunt illuminated by the focused 
laser beam. Voltages were read every 0.1 mm in the 
x-direction which provided adequate resolution. Upon 
reaching the end of the sample, the x-position was re- 



turned to zero and the y-position was incremented by 
0.1 mm. The program Labview was used to control the 
sample position and to acquire the data. The experi- 
ments were performed under the same conditions for the 
bare sample and the MSH in order to produce consistent 
results. The position-dependent data acquisition, as out- 
lined above, was repeated for various temperatures, bias 
currents, and power densities. 



III. THEORETICAL MODEL 

When the laser illuminates a circular region of radius 
Tg, at a position {xi,yi), electron-hole pairs are created 
within the region and the individual carriers diffuse away 
from the spot center. As noted by Dember-''^ many years 
ago, the electrons, which have a much higher mobility 
than the holes {ji^/jih ^ 20 in GaAs)fi& diffuse away 
more rapidly from the focal spot region which causes 
an excess steady state positive charge to develop there. 
Outside the focal region, the diffusing carriers eventually 
recombine due to a number of processes. Nevertheless, 
the excess charge in the focal region perturbs the poten- 
tial at the voltage probes and this perturbation depends, 
not only on the position of the focal spot but also on a 
number of other factors, such as temperature, excitation 
wavelength, power density, etc. Most importantly, the 
voltage probe perturbation also depends on the proxim- 
ity to the shunt that, if present, enhances the removal of 
the photo- induced carriers and thus quantitatively affects 
the magnitude of the excess charge in the focal region. 
The processes described above can be modeled using the 
semiconductor drift-diffusion equations with terms that 
account for photo-generation and recombination. The 
electric field appearing in these equations is obtained 
from the solution of the Poissons equation with the excess 
carriers as the source terms. Because the skin depth (100 
nm at a wavelength of 476.5 nm) and the hole diffusion 
length computed according to Chuang^^ [Lh = 0.1 nm 
(see also discussion in section III. A. below)] that cor- 
responds to several carrier lifetimes is a fraction of our 
sample thickness, the problem becomes two dimensional. 



A. The Drift-Diffusion Equations 

Following the arguments of McKelvej*ii and Chuang,— 
when photo-generated carriers are present, the current 
has the usual resistive and diffusive terms, 

J„ = qifinnE + DnVn) (2) 
Jp = qif^pPE - DpWp), (3) 

where n and p represent the electron and hole number 
densities and fin and fip are their respective mobilities. 
The quantity q = \e\ is the magnitude of the electron 
charge. The quantities D„ and Dp are the diffusion co- 
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efHcients defined by the Einstein relations, 



q 



(4) 
(5) 



Tlie effects of recombination and e-li pair generation must 
be taken into account in the equations of continuity so 
that we have 

-V-Jn + G{x,y)-Rn{x,y) = ^ (6) 

--V-Jp + G(x,v)~Rp{x,y) = ^, (7) 

where G{x,y) is the spatially dependent generation rate 
per unit volume, and Rn{x, y), Rp{x, y) are the recombi- 
nation rates for electrons and holes. 

Using Eqs. (2) and (3) in the above equations gives 
the steady state semiconductor drift-diffusion equations, 

V- {Dr,Vn + finuE) + G{x,y) ^ Rnix,y) = (8) 
^■{Dp^p- fLppE) + G{x,y)-Rp{x,y) - (9) 

The carrier number densities can be written as n = rio+n 
and p = po + p. In the case at hand, the materials are 
n-type with po = 0, and the equilibrium electron number 
density rio is much larger than the excess electron density 
h. As noted above, in GaAs the electron mobility is 
roughly 20 times larger than that of holes. As a result, 
the excess electrons (n) diffuse away rapidly from the 
laser spot and therefore can be viewed as a small increase 
in the large equilibrium electron density. Moreover, since 
Po = 0, the hole number density is simply equal to the 
excess hole number density p. From these assumptions 
follows an approximate charge neutrality, in which the 
charge density remains close to uq. The drift-diffusion 
equation for holes can then be written as, 

W-{DpWp- HppE) + G{x,y)-Rp{x,y) = 0, (10) 

where the quantity E is the local electric field. The re- 
combination rate for holes is assumed to be proportional 
to the hole number density, Rp = SilEjUl^ where r = 0.5 
ns is the carrier recombination timc^^ used in our cal- 
culations. The generation rate G{x,y) is assumed to be 
a Gaussian, with the average radius being equal to the 
spot radius of the illuminating laser 



G{x,y) = GoBxp 



(11) 



Here Go is the generation rate per unit volume at the 
spot center, and the value of the coordinate r is relative 
to the spot center position {xt,yi). The drift diffusion 
equation for holes is now given by 



V- {DpS/p — fMppE) + GoSxp 



4V./ 



-=0, (12) 
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FIG. 3: (Color online) Excess hole density in (a) the bare 
device and (b) the MSH device at the laser spot position 
(xi,yi) = (4nim, 0.7nim). 



In order to solve Eq. (12), we set (a) p = on the 
current ports 1 and 4, and (b) p — Qvn the metallic shunt 
region. The first boundary condition asserts that any net 
hole distribution on the current ports will recombine with 
the electrons in the Indium leads. The second boundary 
condition is due to the fact that the shunt is metallic 
and therefore any hole states will be short lived. These 
boundary conditions are straightforward to implement 
within the framework of the FEM. Figure 3 shows the 
hole distribution for both the bare and shunted device. 
The hole distribution for the bare device is delocalized as 
compared to that of the shunted device. 



B. The Finite Element Method 

The solution of Eq. (12) is obtained using the Galerkin 
FEMfi^iSS with the weight functions being the same func- 
tions as the interpolation polynomials used to represent 
the number density p in each finite element. Multiply- 
ing Eq. (12) by the Galerkin weight function 0(a;, y) and 
integrating gives 



(jN- {DpVp ~ ^ippE)dA - 
-Go II (pexp 



n / r 
4 



dA. (13) 
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The first integral in the above equation is integrated by 
parts to yield 



- V^- {DpVp - HppE)dA + 
(j){DpS/p — upE)- hdl 

'o jj <l>exp 



TT 1 r 

'4 



^dA = 



dA. 



The second term, the surface term, in the above equation 
represents a flux of holes through the external boundary 
and is included only at the current ports in Fig. 1(c). 
The hole current at the current ports can be calculated 
using 



^port 



qd J DpVp- hdl, 

port 



(14) 



where d is the thickness of the sample. These diffusive 

currents must be included in determining the net current 
flowing through the device. The applied voltage neces- 
sary for maintaining the net current at its given value 
is calculated by determining the effective device resis- 
tance including all carriers. This voltage is applied as 
the boundary condition for the Poisson equation, 



V- (eVV) 



-QP, 



(15) 



where e is the material permittivity. Equation (15) is 
solved using the Galerkin FEM for the electrostatic po- 
tential due to the excess hole density with the boundary 
conditions that the applied voltage mentioned above is 
applied across the current ports 1 and 4, while the volt- 
ages at ports 2 and 3 (see Fig. 1) are determined by the 
solution. 

The steps for calculating the potential due to the 
photo-generated carriers are outlined below for one laser 
spot position: 

1. Solve the drift-diffusion equation for p = (Eq. 
12). 

2. Calculate the current on the input and output cur- 
rent ports (Eq. 14). 

3. Determine the applied voltage V14 from the require- 
ment that the net current be the given current 
through the device. 

4. Solve Eq. (15) with V14 applied across the current 
ports. 

5. Calculate the voltage difference V23 between the 
voltage ports using the solution of Eq. (15) from 
the last step. 



IV. RESULTS AND DISCUSSION 

A. Dependence of EOC on Laser Spot Position 

The measured voltage V23 and theoretical results (solid 
lines) are shown in Fig. 4 for y( = (0.2 mm and 0.7 mm) 
as a function of X£ for both the hybrid and bare samples 
at a temperature of 15 K with a power density of 6.3 x 10** 
W/cm^. The theoretical shape of the hole distribution is 
governed by the mobility, and thus the mean free path. 
The mean free path is calculated from 



I = VpT 



= /i\^^(^), (16) 



which takes into account the Fermi velocity, vp, Fermi 
wave vector, kp, scattering time, r, effective mass, m*, 
mobility, fi, thickness, t, of the sample, and the vol- 
umetric carrier density, n. As increases there is a 
peak (valley) corresponding to the positive V2 (negative 
V3) voltage lead. As is evident in Fig. 4, the theoret- 
ical model reasonably approximates the positional de- 
pendence along the x-direction. The generation rate Go, 
which was the only adjustable parameter, was varied such 
that the measured and calculated voltages coincided at 
the peak value. The value of the voltage at other .T£ po- 
sitions was determined using the same generation rate 
which was used to match the peak value. This process of 
choosing the generation rate was done separately for each 
value of yi and for both the bare and shunted sample. 

Figure 5 shows the voltage V23 for a;^ = 3.3 mm as a 
function of yi, for both the hybrid and bare sample with 
theory represented by solid and dashed lines. As yi in- 
creases, the absolute value of the voltage approaches zero 
for both the hybrid and bare samples. In addition, with 
increasing yi, the charge distributions are further away 
from the voltage leads and closer to the shunt, thereby 
lowering the measured voltage. As a result, the voltage 
in the hybrid sample decreases as the laser spot is moved 
closer to the shunt. Because of these effects, both V23 and 
the EOC decrease as a function of the y position, with its 
peak being at yi = 0. Similar to the xi dependence, one 
generation rate for the ye dependence was chosen such 
that the theoretical voltages matched the experimental 
values over the entire range of ye values. 



B. EOC Dependence on the Bias Current 

The voltages in both the hybrid and bare samples show 
an offset such that when the laser is not illuminating the 
sample, the measured voltage is nonzero. The voltage 
offset was seen to be additive to the voltage produced by 
the excess carriers over the range of bias current studied. 
In order to calculate the EOC, the offset was removed 
for both the bare sample and the MSH so that the data 
reflects the effect of the excess carriers alone. In addition, 
the voltage offset was found to be proportional to the 
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FIG. 4: (Color online) Voltage as the laser spot is moved 

in the x-direction for ye = 0.2 mm (Q, ) and ye — 0.7 

mm (□, — ) the bare sample (a) and the MSHs (b). The curves 
represent theoretical calculations and the points represent ex- 
perimental measurements. Experimental measurements were 
taken at 15 K with a power density of 6.3 x 10^ W/cm^. 
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FIG. 5: (Color online) Voltage as the laser spot is moved in 
the y-direction for ye — 3.3 mm for the bare sample (□, — ) 
and MSHs (Q, ). The curves represent theoretical calcu- 
lations and the points represent experimental measurements. 
Relevant experimental parameters are given in Fig. 4 caption. 



nation of carriers. Note that the range of power spans 
approximately two orders of magnitude. 



D. EOC and Temperature Dependence 



bias current which suggests that it is associated with the 
intrinsic sample resistance and not the perturbation that 
we are interested in. 



C. EOC and Laser Power Dependence 

The graph in Fig. 6 shows a plot of the dependence 
of V23 on the illuminating optical power density P. The 
theoretical curves indicate that the measured voltage de- 
pends linearly on the generation rate. Similarly, the ex- 
perimental data indicates that the measured voltage de- 
pends linearly on the optical power density. These two 
observations suggest that the generation rate is propor- 
tional to the incident power, as should be expected in this 
regime. For the bare device, the coefficient of proportion- 
ality is atare = 4 X 10^^ J~-^cm~-^, and for the shunted 
device the coefficient is ashunt = 1 x 10^^ J^^cm"^. The 
difference in a for the two devices arises from geometric 
contributions which affect the generation and recombi- 



As reported^ previously, the EOC reaches a maximum 
of almost 500% at 30 K for optimal {xi.yi) values. The 
EOC also displays an inverse relationship with temper- 
ature. Here we will quantify and further elucidate that 
relationship. 

Because the GaAs is degenerately doped, the equi- 
librium carrier concentrations {no,Po) change minimally 
over the temperature range studied.^- The electron 
mobility^^ and dielectric constant are also assumed 
constant over that temperature range. Implicit in the 
diffusion and resistive terms in Eq. 3 is a temperature 
dependence. According to Lovcjoy et alm^ the hole mo- 
bility, /ip, has an inverse temperature dependence that 
can be fit as fip — xT~^ where x is a fixed parameter 
used to fit the data. In Eq. 5, we see an explicit lin- 
ear temperature dependence of the diffusion coefficient. 
Thus the combined temperature dependence from Eq. 5 
yields Dp ~ T^a. As a result, the effect of the diffusive 
component decreases with increasing temperature. 

Figure 7 shows the temperature dependence of the 
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FIG. 6: (Color online) Voltage versus laser power at spot 
position {yi,X() — (3.3mm, 0mm) for the bare sample (□, — 
) and MSHs (Q, - - -)• The lines represent theoretical 
calculations and the points represent experimental measure- 
ments. The voltage has a linear dependence on the laser 
power and the generation rate. The coefficient of propor- 
tionality is abare = 4 X fO^* J~^cm~^ for the bare device and 
Oshunt = 1 X 10^^ J^^cm^^ for the shunted device. Experi- 
mental measurements were taken at 300 K. 
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FIG. 7: (Color online) The temperature dependence of the 
EOC as defined by Eq. (1). The measurements were acquired 
at spot positions xi = 3.3 mm with yi = mm, {{}, ■ ■ •); 

yi = 0.2 mm, (+, ■ ■ ■); yi = 0.4 mm, (□, ); ye = 0.6 mm, 

(x,solid); yi = 0.8 mm, (V, ) ye = 1.0 mm, (*,- •■ -); 

and ye = 1.2 mm, (Q, ^ ^ ^) Experimental measurements 
were taken with a power density of 6.3 x 10'' W/cm^. 



E. EOC and the lateral photovoltaic effect 



EOC for various yi values with a power density of 
6.3 x 10'' W/cm'^ at A = 476.5 nm and xg = 3.3 mm 
with the points obtained from experiment and the curves 
calculated from the model. In the model, the genera- 
tion rate was held constant for the bare GaAs sample, 
while Go for the MSH was scaled to fit the EOC value 
at one temperature. With Go fixed for each yi in this 
manner, the temperature dependence of the EOC was 
then calculated by simply incrementing the temperature. 
This temperature change, as discussed above, was incor- 
porated directly into the diffusion coefficients as well as 
the hole mobility. With the exception of = mm and 
0.2 mm, the EOC was fitted at 15K. These two yg val- 
ues were fitted at 30 K as their 15 K data (shown) was 
deemed errant owing to the difficulty of establishing the 
ye = position, as discussed previously)^ 



Because of the similiarity in experimental setup and 
positional response to laser illumination between EOC 
and the well known lateral photovoltaic effect (LPE),=^ 
a discussion highlighting the differences is appropriate. 
Critical to the discussion of EOC and LPE is the illu- 
mination of the junction. The lateral photovoltaic effect 
requires illumination of the junction between a (doped) 
semiconductor and a metal or a heavily doped (conduct- 
ing) semiconductor, using a configuration in which the 
illuminating beam is incident normal to the plane of the 
junction. In contrast, in the EOC effect, the junction is 
NOT illuminated by the incident laser beam, which is in- 
cident parallel to the plane of the junction. Also, the LPE 
relies on a Schottky barrier at the metal-semiconductor 
interface, while for the EOC devices shown here, the in- 
terface is an ohmic contact. 

Beginning with WallmarkSi all theoretical analyses of 
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the LPE are based on carrier dynamics derived from the 
photocurrent generated across the junction whereas our 
analysis, as explained above, is based primarily on the 
Dember effect and provides a very good fit to the posi- 
tion dependent voltage (see Figs. 4 & 5). Thus, despite 
the similarity of the positional dependence of the voltage 
response seen in some LPE measurements and ours, it 
should be clear that the EOC is a different phenomenon. 

In the LPE, the output voltage at any lateral posi- 
tion of the illumination INCREASES with decreasing 
temperature. ^''•^'^ In the EOC effect in GaAs-In, the out- 
put voltage DECREASES with decreasing temperature. 
Our theoretical model provides a good fit to the temper- 
ature dependence of the EOC (see Fig. 7). This further 
supports the fact that EOC and LPE are fundamentally 
different position dependent effects. 

V. CONCLUSION 

The drift diffusion model of EOC has been successfully 
employed to fit experimental data taken using GaAs-In 
MSH structures. The model based on the FEM closely 
fits the voltage data and the resulting EOC calculation. 
The position dependence in x and y is incorporated in 
the model via the drift diffusion equation for holes based 
on the mean free path, and thus the mobility. By includ- 
ing a generation and recombination rate a reasonable fit 
was obtained. The model also reproduced the decrease 
in EOC with increasing yi. The EOC decreases with 
increasing temperature T, which the model directly ac- 
counts for through the temperature dependence of the 
mobility and the diffusion coefficients. 

As mentioned briefly in our preliminary report^ EOC 



is similar to but fundamentally different from the well 
known lateral photovoltaic effect (LPE). The main dif- 
ference lies in the fact that the interface is not illumi- 
nated in EOC, whereas in the LPE such illumination is 
critical. Moreover, because of the many factors influenc- 
ing the EOC, it presents a unique opportunity to study 
and model the dynamics of the metal-semiconductor in- 
terface, carrier transport phenomena, and the effect of 
introducing a shunt. 

Because the carrier mobilities are critical in the Dem- 
ber effect, InSb should yield a two orders of magnitude 
increase in the observed EOC, possibly even room tem- 
perature EOC, as compared to GaAs. Using the drift 
diffusion model, we anticipate modeling the EOC of this 
direct gap semiconductor and others. We have not ex- 
perimentally or theoretically optimized for wavelength, 
geometry of the device, or material parameters. With 
the new ability to model EOC as shown, predictions can 
be made to maximize the gain of an EOC device. Once 
optimized, we may expect MSH devices based on EOC to 
impact optical sensing technology, particularly position 
sensitive detectors. 
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